Shear viscosity: velocity gradient as a constraint on wave function 
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By viewing a velocity gradient in a fluid as an internal disturbance and treating it as a constraint 
on the wave function of a system, a linear evolution equation for the wave function is obtained 
from the Lagrange multiplier method. It allows us to define the microscopic response to a velocity 
gradient in a pure state. Taking a spatial coarse-graining average over this microscopic response 
and averaging it over admissible initial states, we achieve the observed macroscopic response and 
transport coefficient. In this scheme, temporal coarse-graining is not needed. The dissipation 
caused by a velocity gradient depends on the square of initial occupation probability, whereas the 
dissipation caused by a mechanical perturbation depends on the initial occupation probability itself. 
We apply the method of variation of constants to solve the time-dependent Schrodinger equation 
with constraints. The various time scales appearing in the momentum transport are estimated. The 
relation between the present work and previous theories is discussed. 

PACS numbers: 05.60.-k, 05.70.Ln, 05.20.Jj. 
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I. INTRODUCTION 

To express macroscopic observables in terms of the mi- 
croscopic parameters of composite particles, one has to 
reduce an infinite number of mechanical degrees of free- 
dom to a few hydrodynamical degrees of freedom[l(. Two 
typical techniques are projection and coarse-graining. In 
the projection method, one can either project the density 
matrix to a reduced density matrix for the relevant de- 
grees of freedom (Schrodinger picture) @, H| or project a 
sum of microscopic mechanical variables to a macroscopic 
slow variable (Heisenberg picture) 0-Q • While in the 
coarse-graining method, excepting the obligatory ensem- 
ble average, one still has to carry out at least one of three 
average procedures: (1) coarse-graining in time @, H|; (2) 
coarse-graining in space [9l-[l2j; or (3) coarse-graining in 
the eigenvalue spectrum of energy and other collective 
variables [1, |1] . 

In the coarse-graining techniques, the spatial coarse 
grained average plays an essential role in formulat- 
ing all macroscopic problems. Macroscopically, a non- 
equilibrium process in a fluid is described by some local 
hydrodynamic (thermodynamic, macroscopic) variables: 
temperature, stress, concentration and macroscopic ve- 
locity field etc[l|. The fundamental assumption behind 
this macroscopic description is that there exists a local 
equilibrium at each coarse-grained spatial 'point' and at 
every coarse-grained temporal 'moment' [l]]. Similarly in 
the kinetic theories, at least three of the arguments of the 
single-particle distribution function, time, position and 
translational momentum are defined in a coarse-grained 
sense (far more coarse-grained than the requirement of 
uncertainty principle) [J . As an alternative to the 
distribution function in phase space, Onsager introduced 
a distribution function for the coarse grained moments of 
the conserved quantities (displacements) [H], EH ■ To ob- 
tain hydrodynamic equations and transport coefficients, 
M. S. Green noticed that the resolution of a macro- 
scopic measurement is finite, so that one must eliminate 
the high wave number components of a collective vari- 



able which are are unmeasurable. The observed value 
of a hydrodynamic variable is given by averaging the 
truncated [T3, EJ collective variable over the distribution 
function of Markoff process [lril [l7j. 

It is well-known that temporal coarse-graining (TCG) 
is not necessary for deriving the macroscopic Maxwell 
equations (MAME) from the microscopic ones[9l-[TT|. A 
safe lower limit of length Lq for a macroscopic descrip- 
tion of reflection and refraction of visible light can be 
taken as I0 2 A 12]. The time scale to associated with 
Lq is in the range of atomic or molecular motions, a 
TCG at a time scale to is meaningless [l2| . Moreover, 
if one averages the motion of the particles over a time 
scale longer than to, the scattering phenomena will be 
smeared [lj|. In the MAME, the sources of fields are the 
total charge density and the total current densityfT^lrsj. 
At the same time, the induced charge density and the in- 
duced current density are also the responses of the system 
to the external fields. Since no TCG is taken in deriv- 
ing MAME, one should not require TCG in computing 
current density (conductivity) [19|-|22| . The energy dissi- 
pation of system is taken care of by the thermal contact 
with a thermal bath ^j^j . It seems that TCG is not 
needed for the irreversible processes caused by the me- 
chanical distrubancesfl-El, [HM HI, HI . 

An open question is whether TCG is necessary to 
describe an irreversible process induced by an internal 
(thermal) disturbance[25|. Several schemes have been 
designed to replace an internal disturbance with a fic- 
titious mechanical disturbance (27l - [3l| . But these meth- 
ods assume a priori that the Navier-Stokes equation is 
valid[32[. In other theories, TCG is often taken along 
with s pati al coarse-graining and ensemble average H 0, 
SHElHlll. In particularly, TCG is viewed as a criti- 
cal step producing irreversibility in the kinetic approach 
(master or Boltzmann equation) 0, [24[ . 

Although a velocity gradient is often viewed as an in- 
ternal disturbance, it can be realized in a 'mechanical' 
manner. In a fluid, two obvious ways to produce a ve- 
locity gradient are (i) moving a boundary plate which 
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confines the fluid; and (ii) stirring the fluid with a rod 
in the middle. In both situations, a velocity gradi- 
ent is produced by the interaction between the fluid 
and an external object at the solid-fluid interface. A 
non-equilibrium 'Hamiltonian' has been phenomenologi- 
cally introduced to effectively compute viscosity [33l-l35| . 
In addition, for a mechanical disturbance [1^, |26| . one 
may define microscopic response in a pure state. The 
observed macroscopic response (consequently transport 
coefficient) is given by spatial coarse graining the mi- 
croscopic resp onse and averaging over possible initial 
conditions [U HO, [23|, |36j]. It is worthwhile to explore 
whether we could describe viscosity in a more 'mechan- 
ical' manner: (i) not invoking temporal coarse-graining; 
and (ii) not assuming local equilibrium as the starting 
point. 

In this paper, we show that shear viscosity can be 
strictly calculated from a time-dependent Schrodinger 
equation, and time coarse-graining can be avoided. One 
may repeat the same procedure for bulk viscosity mu- 
tatis mutandis. In the discussion, we will only use the 
Schrodinger picture. Because a hydrodynamic quantity 
is a bilinear form of the many-body wave function, all 
the conclusions are valid for any set of identical parti- 
cles: bosons or fermions. Only when we take concrete 
approximations for the many-body wave functions, do 
we need to know whether the particles are fermions or 
bosons. 

In Sec. HH we first express the velocity gradient in a 
fluid as a constraint on the wave function of the fluid, cf. 
Eq. (flB")) . To prescribe an internal disturbance like veloc- 
ity gradient, one needs the wave function of the system. 
The wave function must be determined self-consistently 
with the given velocity gradient. With velocity gradient 
as a constraint, the evolution equation for wave function 
is then derived using the Lagrange multiplier method[37|. 
cf. Eqs. (l23l24l25[) . From the appearance, it is a time- 
dependent Schrodinger equation. Two additional terms 
[Eqs. ([241251) ] appear in the 'Hamiltonian'. Both of them 
contain the occupation probabilities of admissible initial 
states, which reflects the fact that a velocity gradient is 
an internal disturbance. On the other hand a mechanical 
perturbation is completely specified by the time depen- 
dence of external field 26]. 

In Sec. IIII1 we first discuss the entropy production rate 
of the system + bath, cf. Eqs. (|33l32l) . Secondly, we check 
the mass conservation law, cf. Eqs. (|34l35p . Thirdly, the 
microscopic response Eq. ([37|) to the velocity gradient (the 
momentum flux in a pure state) is obtained from the mo- 
mentum conservation law in a pure state. The dissipa- 
tion caused by a mechanical perturbation depends on the 
occupation probabilities of admissible initial states. Be- 
cause a velocity gradient is an internal disturbance, the 
dissipation caused by a velocity gradient depends on the 
squares of the occupation probabilities of admissible ini- 
tial states. This feature is clearly seen in the time rates 
of change of energy, mass and momentum. 

Representing velocity gradient as a constraint on the 



wave function [Eq. (fT6l) ] depends critically on the assump- 
tion that a suitable spatial coarse-graining is adequate 0- 
[l2| to describe internal friction. In Sec. IIVI we show that 
the spatial coarse-graining average automatically con- 
tains a coarse-graining in time [Eq.([42"j)] and a coarse- 
graining in eigenvalue spectrum of collective variables 
[Sec. uTB] . 

In Sec. El we apply the method of variation con- 
stants to obtain the Lagrange multipliers which appear 
in the solution of the time-dependent Schrodinger equa- 
tion. The Lagrange multipliers characterize the interac- 
tion on the system exerted by the boundary plates. The 
phenomenological non-equilibrium 'Hamiltonian' (33l - [35| 
can be obtained from Hd, the interaction of system with 
plates linear in velocity gradient. The shear viscosity is 
read from the macroscopic momentum flux (|59p which is 
deduced from the wave function of system at some mo- 
ment. In Sec. I VII we explain that applying Dirac pertur- 
bation theory to transport process is valid. Various time 
scales in the momentum transport process are discussed. 
We show that the length scale of spatial coarse-graining 
is determined by an intrinsic time scale. 

In Sec. IVII1 we show that the macroscopic stress ten- 
sor derived from the microscopic response method is the 
same as that derived from averaging the momentum flux 
operator over the density matrix. We developed a cu- 
mulant expansion for the density operator. When we re- 
place the operator in the exponent with its expectation 
value, we reproduce the non-equilibrium density matrix 
obtained from other approaches. 

II. EVOLUTION OF STATE DRIVEN BY 
VELOCITY GRADIENT 

A mechanical disturbance on a system can be ex- 
pressed by the coupling between the mechanical degrees 
of freedom of the system and some specified external pa- 
rameters which may depend on time and position. The 
disturbance caused by a velocity gradient is more com- 
plicated. We will see that to describe a fluid with a ve- 
locity gradient, the occupation probabilities of admissi- 
ble initial states enter the 'Hamiltonian' of the fluid. In 
Sec lIIAl we briefly summarize the microscopic response 
method for a mechanical perturbation [l9l |20L |23| . [36| . 
Later, we will show that one can still define a momen- 
tum flux in a pure state (microscopic response to veloc- 
ity gradient), viscosity can be calculated from a modified 
time-dependent Schrodinger equation. 

A. Macroscopic response to a mechanical 
disturbance 

Consider a N— particle system S described by Hamil- 
tonian Hq, we use indexes a,b, ■ ■ ■ to label the eigenvalues 
{E a } and eigenstates {3> a } of Hq: 

H <S> a = E a <P a . (1) 
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If the system S is in thermal contact with a thermal 
reservoir B at temperature T, and is in equilibrium with 
it [13, HH, the probability that the system is in state $6 
is 

P * 6 = e -E b /k BT/ ^ e -W_ (2 ) 

a 

Because the system is macroscopic, comparing the en- 
ergy of system S, the energy exchange rate between the 
system and the bath is slow{38j. In a time period much 
longer than a macroscopic measurement, the system can 
be viewed as isolated, and described by a wave function. 
For a macroscopic fluid, the detailed dynamics of the sur- 
rounding bath is not important. 

If a mechanical disturbance described by Hamiltonian 
Hf is exerted on S, the state "J of S will evolve according 
to 

ihdV/dt=(H +H f )y. (3) 

The above description requires that [23|, [24J S is in good 
thermal contact with B such that the energy generated 
by the disturbance Hf inside S can be rapidly transferred 
into £?, and the system is maintained at the temperature 
T of bath B. If we assume that the system is initially in 
state $ a , the state ^ a (t) of system at time t is determined 
by Eq.Q with initial condition ty a (t = — oo) = $ a . For 
convenience we adiabatically introduced Hf. By means 
of the microscopic conservation law in a pure state, we 
can define^ the corresponding microscopic response 
(flux) in state \& a (t). 

One can always spatially resolve a quantity A^ a (R, t) 
in state \l/ a into its Fourier components: 

Because any instrument has finite spatial resolution, one 
can only detect those A^ a (k, t) with |k| < ho, where kg i s 
some instrument limited cut-off wave number^. [Til [l7j. 
To express a macroscopic measurement, one has to make 
spatial coarse-grained average fill Il2| : 

A C *»M) = J dR/(r-R)yl*»(R,t), (5) 

where /(s) is any function satisfying: (1) J <fs/(s) = 1; 
and (2) the Fourier components of /(s) tends rapidly to 
zero for |k| > kg. The integral in Eq.© is over the 
sample. 

Since the system is in thermal contact with a reservoir, 
i.e. a member of a canonical ensemble, one does not 
know which state the system is initially in. The measured 
quantity is an ensemble average of A* aC (r,i): 

A{v,t) = Y,P* a A c **{v,t), (6) 

where P* a (— oo) = _P* a is the probability that the sys- 
tem is in state ty a (t = — oo) = $ a in the remote past. 



Since the microscopic response A^" (R, t) is a bilinear 
form of it can be computed from Eq.@ rather than 
from the Liouville equation for the density matrix. Since 
the spatial coarse-graining (|5|) and the average over ini- 
tial conditions (|6]) do not involve the temporal evolu- 
tion of the system, the macroscopic response A(r, t) and 
consequently transport coefficient are also determine d by 
the time-dependent Schrodinger eauation[19l l20l. \24. 136| . 
The averaging procedure (|5I6[) is simpler than, and equiv- 
alent to Kubo's linear response theory [36} . 

The procedure (|5I6[) aims at the steady states {^ a {t)} 
of S under a monochromatic driving. The general lin- 
ear casual response to a disturbance with arbitrary time 
dependence will be discussed in Sec lVIB2l 

B. Velocity gradient as a constraint on wave 
function 

Consider a fluid composed of N identical particles with 
mass m. We use indexes j, k, I ■ ■ ■ to label the particles. 
Denote the interaction energy between two particles at 
rfc and r; as f/(rfe,r;). Then the Hamiltonian of system 
is 

k=l k,l(kjil) 

Put the system in an external field, denote the interac- 
tion potential energy for a particle at in the external 
field as V(rk,t). The interaction Hamiltonian with ex- 
ternal field is 

N 

H f = Y, v ^- ( g ) 

fe=l 

In the remote past t = — oo, suppose that the system 
is in state $ a - Driving the system with a velocity gradi- 
ent and Hf, $ a will evolve into a state ^(rx • • ■ rjy; t) at 
time t. Eq.([2]) indicates that we do not have the full 
knowledge of the system at the initial moment. The 
evolution of the system can then be understood to be 
probabilistic [1^]. The interaction U among molecules 
renders the macroscopic motion involving a velocity gra- 
dient irreversible [39j. 

1. Macroscopic velocity field in a fluid 

Using the definition of mass density operator || [I?} 
p(r) = Ylj=i m5(r — rj), the mass density for state ^ a (t) 
is a function of time and position: 

p**(T,t) = Nm[dn*Z* a , (9) 

where dr\ = d,Y2 ■ ■ ■ dr jv , the arguments of ^ a are 
(r, T2, f3, • • • , tat; t). The momentum density operator is 
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FIG. 1: Experimental setup for velocity gradient: the top 
plates ABCD moves toward +x direction (right) with speed 
v x0 . The bottom OEFG, front (right) bank AEFD and back 
(left) bank BOGC are static. The width of channel is l z , the 
depth of channel is l y . 



denned as0,[T3| 




(10) 

where 7 = x, y, z are indexes for three Cartesian compo- 
nents. The momentum density in state a is 

According to Eqs. (|5I6[) . the measured macroscopic 
mass density is 

p(r,t) = ^P*» fdRf(r-R)Nm f dr^ a %, (12) 

where the arguments of ^f a are (R, r 2 , r 3 , ■ • • , rjv; t). The 
measured macroscopic momentum density is 

S 7 (r,t) = / dR/(r-R)S*«(R,i). (13) 

The macroscopic velocity field should be defined asffill7l. 

vy(r,t) = Sj(r,t)/p(r,t), j = x,y,z. (14) 

It satisfies the requirement that the macroscopic velocity 
is the macroscopic momentum of unit mass of fluid 40]. 

2. Velocity gradient 

We first describe the simplest experimental setup for 
measuring shear viscosity in Figjl] Let the fluid flow in 



a rectangular channel along the x direction. We choose 
a coordinate system such that the two banks BOGC and 
AEFD of the channel are z = and z = l z , the bottom 
surface OEFG of channel is y = 0, the top surface ABCD 
of the fluid is bordered by a plate y = l y with velocity 
v x q. The points on the bottom plate OEFG are denoted 
as r*- 1 ' (x, 0, z) with — 00 < x < 00 and < z < l z , 
the boundary condition on bottom is u 7 (rW,t) = 0, 
7 = x, y, z. The points on the top plate ABCD are de- 
noted as r( 2 \x, l y , z) with —00 < x < 00 and < z < l z , 
the boundary condition on the top plate is v x 
v x o, v y{ r 1 t) = v z {r^ 2 \t) — 0. The points on the back 
(left) bank BOGC of channel are denoted as r^(x, y, 0) 
with — 00 < x < 00 and < y < l y , the boundary 
condition on the back bank BOGC is u 7 (r^,t) = 0, 
7 = x, y, z. The points on the front (right) bank AEFD 
of channel are denoted as r^ 4 \x, y, l z ) with —00 < x < 00 
and < y < l y , the boundary condition on the front bank 
AEFD is u 7 (r( 4 ),i) = 0, 7 = x, y, z. These no-slip condi- 
tions for a viscous fluid are the boundary conditions for 
the Navier-Stokes equations [4Qj. We imposed a velocity 
gradient v x o/l y in the fluid. 



3. Velocity gradient as constraints on wave function 

To enforce a velocity gradient as a constraint on the 
wave function, it is helpful to allow an arbitrary ve- 
locity field on each boundary surface. We use a two 
dimensional mesh to characterize the positions of the 
points on a boundary surface. For a point on the a th 
plate (a — 1,2,3,4), we use two indexes \i a and v a 
to describe its position, where \i a — 1,2, ••• , Mi CT and 
v„ = 1, 2, • • • , Mia- Thus the position vector of a point 
on the a th plate with indexes and v a is denoted as 
r^lv , there are M\ a x M 2a points on the a th boundary 
surface. A general velocity field on the a th boundary sur- 
face may be specified n 7 (r^l„,t), 7 = x,y,z. In terms 
of Eq. (fl"4"j) , the boundary condition on the a th boundary 
surface is expressed as 

Sy (r£l„ , t) - p(r%> Va , t)v 7 (rW o , t) = 0, (15) 

where v 1 {v^J l/ljl t) is a know function of (Tj^J u<r ,t). Mak- 
ing use of Eqs. (|12lllll3|) , Eq.(Tr5]) becomes 

£ f dR/(r^ a - R)N [ dn (16) 



Eq. ([T6|) represents M c = 3 x 4 x M\ a x M20 constraints 
on the wave function of system. The evolution equation 
of ^ a {t) has to be modified with constraints (JTgJ) . 
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C. Evolution of state in a velocity gradient 

To consider the effect of constraints (fT6|) on the state 
ty a of the many-particle system, we first notice that the 
Schrodinger equation 

iM/ at = (H + fl»*„ (17) 

is the Euler equation SS/S^Z = of the action (37l. lilf: 

S = J'dtj K,V* at ,y* ar .), (18) 

where tf at = d& a /dt, = V rj ¥ , dr = 

dv\dv2 • ■ ■ drpf is the volume element in configurational 
space, the arguments of ^ a are (ri,r 2 , • ■ • , rw;f). The 
well-known Lagrangian density for Eq. ([T7j) ispTjj : 

N h 2 

C = im* a * at - y £— Kr^ar, (19) 
3=1 



where the Lagrange multiplier Xp°v a (t) is a real func- 
tion of time[42| with dimension length 4 -time" 1 . Roughly 

speaking, XjJfX (*) i s a product of the rf h component of 
velocity at point vfc} Va of the a th plate and the disturbed 
volume of fluid by a molecule at TpJ Va . We will see this 
in more detail in Sec IV Al 

The equation of motion for ^f a (t) can be obtained from 
the Euler equation: 

c^l_Vv dC c d dC c 

d** a £j r '5** r . dtd** at - ■ { ' 

It is 

iMat = {Ho + Hf + H e + H d )* a) (23) 

where 

H C = P^ X t:l® /rfR/(r (ff) -R) (24) 



k=l k,l(k^l) 

We apply the Lagrange multiplier methodfllj to find 
the evolution equation of ^ a (t) driven by a velocity gra- 
dient specified by Eq. (|16[) . First of all we symmetrize 
Eq. ([T6"|) respect to all particles. Secondly, we multiply 
the expression with X]JfJ a (t) and integrate it over the 
time interval p2l| [tijtaj: 




-m« 7 (r^,t)* 6 *J]=0, (20) 

where the 6 function is produced when we change J dr\ 
in Eq.(fT6f to / dr in Eq. ([20|) . The macroscopic velocity 
at every point on each plate is a constraint. The overall 
Lagrangian density including the constraint (|20[) is 

N 

/( r ^- R )E*( r i- R ) 

3=1 



N 
3=1 

and 

{ *( r . _ R>fi— + - V r . 4 - (25) 

One may notice that Eq. ([23| contains only ^ a and 
its partial derivatives. Other do not appear. In 
this sense, the velocity gradient is a special mechan- 
ical perturbation. On the other hand, the solution 
* Q (t) of Eq.pg]) involves M c parameters {\]$ B (t)}, 
that should also satisfy M c constraints (|16[) which mix 
all {^(i)}. Thus velocity gradient is also an internal 
disturbance [1, [H[ . 

In H c and Hd, the integrals behind the summa- 
tion signs are the macroscopic momentum density 
and negative momentum density operator respectively. 
Thus the subtracted momentum density automatically 
happens [11] in H c + Hd , which has to be introduced in 
some previous theories Q- 

Because Hd ^ H^, one cannot find a system with 
wave function W({rj};t) = ^f* ({ij}; — t) such that 
ihd^'/dt = (H + H f + H c + H d )^', so that the pro- 
cess described by Eq. (|2"3"|) is irreversible (39j. 

In classical mechanics, constraints can also be dealt 
with the Lagrange multiplier method [42T(45j . A con- 
straint produces a reaction force on each particle in the 
system. The constraints (1161) enforce the interaction on 
the system from four the plates. 
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Defining the macroscopic velocity field of a fluid by 
Eq.([Tl|) has a profound consequence: the resultant con- 
straint (fT6| on wave function is bilinear. The evolution 
equation (|2"3")l for ^ a (t) is linear and does not include 
other states. The evolution of the state of system driven 
by a velocity gradient is almost determined in the frame- 
work of pure mechanics. If we had defined the macro- 
scopic velocity field of fluid as 

u 7 (r, t) = J2 P" a I dRf(r ~ R)<° (R, t), (26) 



with 



•(R,i) = S*«(R,t)/p*»(Rt), (27) 



we would have two serious difficulties: (i) this definition 
breaks the basic requirement that the macroscopic mass 
flux density S* 7 (r, t) must always be the macroscopic mo- 
mentum of unit volume of fluid[40|; (ii) the resulted tem- 
poral evolution equation for ^ a (t) is nonlinear, and ^ a (t) 
couples with other {'J ;,(£)}. Certainly, assuming the fluid 
is incompressible, the labyrinthine evolution equations 
could be linearized and decoupled. 



D. A comparison of velocity gradient and a 
mechanical disturbance 



For a mechanical perturbation, external fields are some 
specified functions of time. The interaction Hamiltonian 
Hf is expressed by the external fields, the coordinate 
and momentum operators of all particles [ill l23l. l25j . The 
states {^ a (t)} of system are not involved in Hf. On the 
other hand, the macroscopic velocity field of a fluid is de- 
fined in Eq. ([H)) through macroscopic mass density ([12")) 
and macroscopic mass flux (|13l) . The later two quanti- 
ties involve the states {^ a (t)} of system, spatial coarse- 
graining and average over various initial conditions. Be- 
fore specifying the velocity gradient, which involves the 
states {^ a (t)} of system, we cannot determine the evolu- 
tion of the state of system driven by a velocity gradient. 

From the charge conservation law in a pure state, one 
can directly use the time-dependent Schrodinger equa- 
tion (J5J) to derivefljl Ho[ the current density j* a (r,t) in 
state (the microscopic response) from the temporal 
change of charge density /0*"(r, t) in state ^ a {t). One 
does not need the quantal operator of current density [36|. 
although the results are the same from the two points 
of view. In contrast, we cannot use mass conservation 
in a pure state to derive mass flux S 1 *" (r, t) in state 
^ a (t). We have to invoke the definition of momentum 
density (fTUj) . In this sense, the velocity gradient is an 
internal disturbance [1, |25|. The states {^ a (t)} of sys- 
tem driven by velocity gradient have to be determined 
by the equation of motion (|2"31 and the constraints p^|) 
self-consistently. 



III. CONSERVATION LAWS IN A PURE STATE 

In this section, we assume that Xj^fJ^ (t) and {^b(t)} 
are known. The time rate of change of energy, mass den- 
sity and momentum density are derived. These relations 

are formal but exact. We will find AjH^v(t) and {^(t)} 
in SecEl 



A. Time rate of change of total energy 

The temporal or spatial change in local energy density 
is closely related to the temperature gradient. Since S 
is a member of the canonical ensemble at temperature 
T, we cannot talk about the time rate of change of local 
energy density. 

For an isolated system S, the total energy is conserved: 



.E*»(t) 



J dTe itE " /H $* a Ho$ a e~ itEa/h = E a . 



(28) 



Consider a more general state of isolated system S 
with initial value ^(0) = ^ e c e $ e , where {c e } are 
some constants. The time evolution of *(0) is ^(t) = 
~^2 e c e <& e e~ lEct / h . From the orthogonality of eigenfunc- 
tions of Ho, it is easy to check that the average energy 
of <£(£) is conserved: 

£* (0) = J d T V* (0)ff„*(0) - ]T \c e \ 2 E e , (29) 



and 



= J dr^{t)H^{t) = J2 \ce\ 2 E e 



(30) 



Now assume that the system S is in good contact with 
a thermal reservoir B in temperature T . Thus, even if 
we apply a velocity gradient and an external field on the 
system, S is maintained at the same temperature T . Be- 
cause S is macroscopic, in a time period much longer 
than a macroscopic measurement, S can be viewed as 
quasi-closed 38] . Therefore we can use a wave function 
^ a {t) to describe the system. The average energy of S 
in state ^ a {t) is defined as 



£*■(<) = / dr**(t)ff„*«(t)- 



(31) 



Because \£ (t) satisfies Ea.(|2"3"|). one can get the time 
rate of change dE^ a (t) /dt of the energy E^ a it) in pure 
state ^ a (t)- The average energy E(t) of S at time t is 
E (t) = E* a P*<*-E* a (*)- The time rate of energy ex- 
change dE(t)/dt between S and B is determined by 



dE{t) 
dt 



J dRV(R, i)V R - 



(32) 
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V X^J /dR/(r^ -R)^-|-Vr- 

^ ^ j j K n*** > 2 ap 7 



The temperature of 5 is maintained at T by the bath 
£?. All the energy generated in the system is transferred 
to B. The entropy S B (t) of the bath B increases with 
time according to 



dS B {t)/dt = T- l dE(t)/dt. 



(33) 



Substituting Eq.flSU) 



that A 



7(<t) 



(*)/(* 



(°0 



into 
R) 



Eq. 

has 



and noticing 
dimension of ve- 



locity, each term is in the form of temperature -1 
x velocity x momentum x velocity gradient, except the 
first and last term. This is consistent with the macro- 
scopic law of entropy production rate[40(. 



+ E Kit J - V R - 



\7(» 



dR/(r (<T) 



-R)N(N-l) 



B. Time rate of change of mass density 



With the help of Eq. ([23|) . one may calculate the 
time rate of change 0, S3 of p^ a (r,t). Further apply- 
ing the averaging procedure (|5l6p to the expression for 
<9j0* a (r, t) j dt, one obtains the macroscopic law of mass 
conservation: 



dp(r,t) 
dt 



E 



dr~ 



0. 



where 



S;(r,t) = S 7 (r,t)-y"dR7(r-R') E A m1 



(34) 



(35) 



Or 



2~, 



where the arguments of \I/ a are (R, r2, •■• , r/v;i). The 
last term is the power of the internal force near the 
boundary surfaces. In classical mechanics, for a group of 
particles, the internal forces do work on the svstem|46|. 
Only for a rigid body, is the work done by internal forces 
zero. It is interesting to notice that only for the regions 
close to the boundary plates [f(r^Ju, ~ R) decays very 



» 



R| > fc ], the internal forces do 



rapidly when 
work. 

In Eq. (|32[) , we can see an obvious distinction between 
a velocity gradient and a mechanical perturbation. The 
first term contains P $a , is the power due to the external 
field. The remaining terms contains [P* a ] 2 , the square 
of the occupation probability of admissible initial states. 
They are the contributions to power due to the friction at 
boundary surfaces. The first P*° comes from prescrib- 
ing the macroscopic boundary condition as constraints on 
the wave function of a pure state. The interaction with 
plates are expressed by H c + Hd in the 'Hamiltonian' of 
system, cf. Eqs. (|24l25l) . Since dE^ a (t)/dt is determined 
by * oi (t), from Eq.flU]), dE^*{t)/dt includes a P*». To 
obtain dE(t) /dt from dE^ a (t) /dt, the average about pos- 
sible initial states (0 introduced the second P* a . 



xf(r (a) 



R')E[ p * 



]V"(R',t). 



In the R.H.S of Ea.([3"5]l, the first term S 7 (r,t) is the 
usual macroscopic mass flux, and is given by Eq. (|13[) . 
The second term is an additional mass flux caused by the 
given velocity gradient. Two truncation functions [111 Il7| 

f(r^Ju a — R') and /(r— R') indicate that the correction is 
nonzero only near the boundary surfaces. Again, [P $a ] 2 
reflects the fact that a velocity gradient is an internal 
disturbance. 



C. Microscopic response: momentum flux in a 
pure state 

Taking the time derivative of Eq. flTTl) . one obtains 
the time rate of change of the momentum density 0, S3 
5* a (r, t) in state ^ a (t). In terms of Eq. ([23| . one has: 

as**.*) +E ggSM = _g^)^ (M)/m 



dt 



dr B 



dr a 



-p*> y az>> ( t ) df ^' v '~ T) 

^— ' ' dr n 



(36) 
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where a, /?, 7 = x, y, z, and the momentum flux 11*^ (r, t) 
in state ^ a is 



(37) 



<9* a d^* a d^l d^ c 



dr a dtp dr a drp 



2 drpdr a 



-N 



N 

J k=1 



In Eq. (j3"7} . the arguments of $ a are (r, r2 ■ • • rjv; t). 
To obtain the divergence form for the second term in 
Eq. (l37l) . we used the fact that the range of interaction 
force is much shorter than[l7l |32| fc^ 1 . 

Although velocity gradient is an internal disturbance, 
by viewing the velocity gradient as a constraint on the 
wave function, one can still define a microscopic re- 
sponse: the momentum flux 11*3 in state ^ a (t). Ap- 
plying Eqs. (|5l6l) . one obtains the measured macroscopic 
momentum flux: 



1 



N 
*=1 



dri 



-]5(Ti-r). (40) 



To obtain the viscosity, one must average operator (|4"0")) 
over a non-equilibrium density matrix (or distribution 
function) p non which accommodates the effects of veloc- 
ity gradient HI, S3, H^-fll- This reminds us of some 
older theories for the response to an electromagnetic field. 
In those theories, the current operator is derived from the 
Schrodinger equation for an isolated system (without ex- 
ternal electromagnetic field). The resulting current den- 
sity missed the vector potential term and broke gauge 

invariancefHHlii]. 

In Eq.(|3T|). * (t) is the solution of Eq.(|23]l rather than 
that of Eq. (|3T)]) . Thus the influence of velocity gradient 
has been included in LT*^ (r, t). If we know ty a (t), Eq. (j3"7} 
directly gives the microscopic response LT*^ (r, f) to a ve- 
locity gradient in state ^(t). We will compare the re- 
sults derived from two method in Sec lVIII At this point, 
we only note that Eq. (f3"T|) reduces to the ordinary con- 
servation law in the bulk, and is also correct near to the 
boundary surfaces. For the effects which are first order 
in velocity gradient, Eq. (|36l37[) gives the same results as 
those based on Eqs. (l39l40p . 



n Q/3 M) = ^p*° 



dR/(r-R)nJ§(R,i). (38) 



After ^ a (t) is obtained, the shear viscosity can be read 
out from Eqs. (|37l38l) . 

In the RHS of Eq. (l36|) . the first term is the momentum 
density produced by the external field per unit time, the 
second term is the momentum density produced by the 
constraint (|16[) per unit time. Because the moving plate 
drives the fluid, S gains momentum and energy. Thus 
there are source terms in the equations for the time rate 
of change, cf. Eqs. (|3"o) and (|3"21) . Since no mass is pro- 
duced by a velocity gradient, there is not a source term 
in the time rate of change of mass density, cf. Eq. (IMl) . 

There is a subtle difference between Eqs. ([361371) and 
the corresponding results in previous theories [171 l32l |47T— 
[54|. In all previous theories, conservation laws and fluxes 
(responses) are derived for an isolated fluid system with- 
out imposing velocity gradient. The system satisfies 

ihd^o/dt = H a ^o (39) 



or Newton's equation for an isolated system 117 
l5(il - [54l ]. The operator of momentum flux is [71.147 



N 1 

iMr) = J2 - r ) - M 



IV. SPATIAL COARSE-GRAINING 

Let us check whether the average procedure (|5l6p is ad- 
equate and enough to describe the dissipation produced 
by a velocity gradient. 



A. Spatial coarse-graining includes temporal 
coarse-graining 

In a liquid, suppose that the force between two neigh- 
boring molecules is U /d, where d is the average distance 
between two molecules, U is the interaction energy be- 
tween two molecules at a distance d. The time needed 
to relax a change h/d in momentum is t m \ ~ K/U. 
The time needed for a molecule to move a distance d is 
t m 2 ~ d/(h/md). Near the equilibrium state, the poten- 
tial energy is balanced by kinetic energy U ~ h 2 /(md 2 ), 
and the two microscopic relaxation times are the same 
order t m i ~ i m 2 and will be denoted as t m . 

We show that the spatial coarse-grained average con- 
tains temporal coarse-graining automatically. Apply- 
ing the spatial coarse-graining average ([5]) to the micro- 
scopic law of mass conservation in a pure state the 
time rate of change of a conserved quantity must be the 
same order as the divergence of the corresponding flux 
dp c *«(r,t)/dt ~ £ 7 dSf«(r,t)/dx 7 , i.e. 



mn h 

— k -n, 

t d 



(41) 
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where n ~ d~ 3 is the number density of fluid. The re- 
laxation time t c i for a fluid drop is t c i ~ md/(k h). 
Because the cut-off wave vector is ko, the time needed 
for a molecule to pass a distance k^ 1 is k^ 1 /(H/md) ~ 
(kod)~ 1 t m 2, the same order as t c \- The temporal coarse- 
graining is automatically realized in the conservation law 
of mass by the spatial coarse-graining ([5]) . 

Applying spatial coarse-graining to the microscopic 
law of momentum conservation in a pure state ^/ a , one 
has dSf'(r,t)/dt ~ Ep9Uf^(t,t)/dr , i.e. 

n h h h 

t d d md ' 

the relaxation time of momentum in a fluid drop is the 
same as t c \. The temporal coarse-graining is automat- 
ically realized in the conservation law of momentum by 
the spatial coarse-graining (0). We do not have energy 
relaxation. 

We derive the coarse-grained time scale t c \ from an- 
other point of view. There are N\ ~ (kod)~ 3 molecules 
in a coarse-grained fluid drop. Because the range of 
molecular interaction is only extended to d, a fluid drop 
interacts with its surrounding only through its surface. 
Thus the force on a fluid drop is N^ 3 U/d. The accelera- 
tion of the fluid drop is a ~ {Nl /3 U/d)/{mN 1 ), the time 

1/3 

needed for the fluid drop to move a distance s — N x ' d 
is t c2 ~ (s/a) 1 / 2 ~ fcf x (m/l7) 1/a ~ (fad) -1 (ft/17). Two 
characteristic times for the coarse-grained fluid drop are 
the same order t c \ ~ t c %, will be denoted as t c . Be- 
cause k^ 1 3> d, we can see t m <C t c , the coarse-grained 
time scale is much longer than the microscopic time scale 
t C 2 ~ {k Q d)~ l t m . 

The time coarse graining induced by the spatial coarse 
graining is fine enough to describe macroscopic mo- 
tion. The shear viscosity is order of rj ~ mU/(hd). 
In a macroscopic measurement, the acceleration on a 
fluid drop produced by the velocity gradient is clm ~ 
^ c|f fco~ 2 [ mn ^cT 3 ]~ 1 - The macroscopic relaxation time is 
order of t M = v x0 /a M ~ Ko(^f d) _1 ]i c ~ tj y /d, where 
l y is a macroscopic length. Therefore tyi 3> t c . In sum- 
mary, we have 

fm<tc<tM- (42) 

For the thermal motion of molecule, we could make a 
similar estimation: the velocity and kinetic energy of a 
molecule should use vt — {k B T/m)^ 2 and k B T respec- 
tively. Relation ([jS)) holds. In conclusion, the spatial 
coarse-graining ([3]) included a temporal coarse-graining, 
which is adequate and enough to describe the dissipation 
produced by a velocity gradient. 

B. Spatial coarse-graining includes coarse-graining 
in the eigenvalue spectrum of the collective variables 

In Sec lII A[ the quantal mean value of any operator in a 
pure state is spatially coarse-grained with Eq.©. Many 



authors choose to coarse-grain the eigenvalue spectrum 
of collective variables. We will illustrate that the spa- 
tial coarse-graining in the microscopic response includes 
coarse-graining in the eigenvalue spectrum of collective 
variables. 

Because the cut-off wave vector is fco, there are iVi ~ 
(kod)~ 3 molecules in a coarse-grained fluid drop. A 
collective variable for a liquid drop is a symmetric 
sum of single-particle variables [HI. The typical inter- 
val A^Vi between two eigenvalues of a collective vari- 
able is Ajvi ~ Aie~ Nl , where Ai is the typical inter- 
val between two eig envalues of the corresponding single- 
particle variable [38|. Excepting macroscopic quantum 
phenomenon (superconductivity and superfluidity), even 
in a microscopic relaxation time t m 2 ~ md 2 /ft, the ac- 
tion for a drop is Ni(h 2 / md 2 ){md 2 / h) ~ (k d)~ 3 h » ft. 
Therefore the collective variables are classical variables. 
Since the relative error for a collective variable is[38| 

— 1/2 

A/// ~ N 1 , the measurement error for the energy of 
a fluid drop is Ajvi ~ Nl /2 {h 2 /md 2 ). The energy interval 
for a single molecules is Ai ~ h 2 /md 2 . A macroscopic 
energy measurement with spatial resolution k$ 1 must in- 
volve N\ e Nl eigenvalues and corresponding eigenvec- 
tors of the energy operator of a liquid drop. 



V. CALCULATING LAGRANGE MULTIPLIERS 

In linear response theory, to obtain the observed 
macroscopic response to a time-dependent mechani- 
cal perturbation, Dirac's method of the variation of 
constants is frequently used to solve the Liouvillc 
equation[26j or Schrodinger eauation [l9l - [2ll [55l |5&!| . We 
first use Dirac's perturbation theory to get Lagrange mul- 
tipliers, and then discuss its applicable conditions. 



A. Lagrange multiplier 

Consider a state ^ a {t) with initial condition in the 
remote past 9 a (t = — oo) = $ a . In first order of H±(t) = 
Hf(t) + H c (t) + Hd-, the normalized wave function is 

tf (t) = e -^$ a (l+ i Im.g( 1 a ))+ J2 sSVK^ 4/ %. 

(43) 

For both j = a and j ^ a, the first order expansion 
coefficient is given by: 

=<&(*)+ E ^:i(tmitha(t). (44) 

The second term in Eq. (|44|) is the contribution from the 
moving plate and three static surfaces: 
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xiV 



ggb 

Si?-, 



0Av bJ 



d<S> b 



dR-y OR 



J dRf(r$ Va - R)N J dn$* c * b , 



(45) 



whcr 



(£y — E a )/h. Denote d as the average 



distance between two neighboring molecules, U) the fre- 
quency of a monochromatic driving velocity v x (r^\t), 



» 



41-1 



characterizes the transition 



caused by a point y)iJv„ at the er plate, and has the 
dimension of time-length -4 . The first term in Eq. (144[) is 
the contribution from the external field: 



dt'e™^ / dT#J[-^^(r* > *')]*a. 



fc=i 



(46) 

The first term in the RHS of Eq. (|4"3l is the free evolution 
term, i1mgaa(t) is caused by the requirement that 1 
j dT^*(t)^ a (t) should be correct to first order in Hi(t). 
In the kinetic theories, one concerns various transition 
probabilities. We only need the j ^ a terms in Eq.(|43|). 

Before we apply the external field and move the top 
plate, the system S is in equilibrium at the bath tem- 
perature T: the system is in various stationary states 
{${,} according to Eq.®. In an ordinary fluid, there 
is no macroscopic mass flux at any boundary surface 

S^(r^,t)=0: 



: E^ b 



R)N 



(47) 



where 



p = m 



+c* db (t)e l{Ed - Eb)t/h W}}, 



J dRf(r^ -R)^P*>N J dn$ b $l (50) 



is the constant mass density of liquid in absent of external 
field and boundary conditions (fl"6|) . 

Making use of Eq. flU} , Eq. (|49| is reduced to a group 
of linear equations for Lagrange multipliers: 

E ^.ZJ%K*)^%(*) = ^(*). ( 51 ) 

r er' er' r a' a' 



where 



J dR/(rW a - R) E P *" N J dT i 



(52) 



bd{b^d) 



,0*1 
■flit, 



) = o, 



where 7 = x, y, 2, er = 1, 2, 3, 4, the arguments of $f, are 
(R, r-2, • ■ • , rjv). Superfluidity is excluded by the require- 
ment (1471) . 

To conserve space in the subsequent calculations, we 
rewrite Eq. (14*31 as 



*«(*) = e-^'^a + c ja (t)e- iE ^ h ^, (48) 



and 



2 v "di^ <9i? 7 



HZ K*) = / - R ) ( 53 ) 



^ . ^ / dn{[^:?]*(t)e- 



where c ja (t) = ilmg£J(t) for j = a; c ja (t) = £fj^(t) for 
j 7^ a. The summation includes both j — a and j a. 
Substitute Eq. (|4"8"|) into boundary conditions (fit*)]), and 
notice Eq.(|4T|) for four boundary surfaces, then we have 

v^r$ Vv ,t)p° =E P * 6 / dRf(r^ Va -R)N J dr x 



(49) 



[ 2 1 d dR 1 dR 7 
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The dimension of B$fJ a is momentum density. The first 



term in Eq. (|52[) is the momentum density produced by a 
point r( CT ) at the a th plate ~ v 1 {r < ^ a \t)md~ 3 . The sec- 
ond term is the momentum density produced by the in- 
terference between the external field and the constraint 
of velocity gradient ~ [(ui d b + w (V d b / d)]d~ 3 , where 
ujf is the frequency of the external field. B a is the net 
momentum density at on the a th plate. The di- 
mension of (A ffIJ )} is momentum density -time -length . 



A 



■j(<t) 7'(ct') 



Lj)d A ] 1 (h/d)d 3 describes the 
effect of a point r^" ' v , at the er'th plate to the point 



>') 



_(«0 



on the a th plate. Eq. ([5"Tj) describes the momentum 
produced by all points on the four 



density at point r^, 



boundary plates. The Lagrange multipliers )C[}fv a are de- 



termined by the coupled linear equations (|5ip . and they 
are in general nonlinear functions of the velocity gradient. 
The lowest order approximation to A^lv is first order 



fJ-cr ' 



,t)/l y : 



,t)/(h/md)][(u; db + Lo)d}d 3 . (54) 



From Eq. ([32j) . we can see that the lowest order dissi- 
pation rate is second order in velocity gradient, which 
is consistent with the result derived from Navier- Stokes 
equation and the second law of thermodynamics [4(|. In 
addition, from Eqs. (|54l24l25l) . H c is second order and 
Hd is first order in velocity gradient. This is similar to 
the coupling between electromagnetic field and material: 
there are both first order and second order terms in the 
vector potential. The small parameter for the pertur- 
bation caused by Hf is (b\V(ri,t)\c)/(H 2 /md 2 ), for the 
perturbation caused by H d is v 1 {r$ Va ,t)/(H/ md) , for the 

perturbation caused by H c is [Vy(r^J Ua ,t) / (h/md)] 2 . 

If we substitute the estimation ([5~4"f into Eq. (l231) . one 
has 



Hd ~ ^{m^frl,! , t){ui d b + u>)d}. 



(55) 



The quantity in {} is the dissipative part of the phe- 
nomenological non-equilibrium Hamiltonian for the j th 
molecule [33|-|35j. Here (o^b + w) corresponds to the strain 
rate (co the frequency of a monochromatic driving), d the 
displacement of an atom, mv 7 {r^} Va , t) the momentum 

of an atom at point rJ^„ CT . Of course, Eq. (j5"5"]) has not 
changed from the boundary surfaces to the bulk[33U35l|. 

From {Xl%\{t)}, Eqs. (|45l46l44p tells us all the expan- 
sion coefficients in Eq. (l4"8l i.e. the state ^ a (0- 



where 



(57) 



L 9i? 7 dRp dR 7 dRp 2 8R B dR. 

N 



-\n J dr^Re-n^-^^SbSl 

is the momentum flux in pure state <I?b- 

By means of Eq. (|48j) . the momentum flux Q37j) in a 
pure state is reduced to 



n^(R,t) = |-iv f dn{J2c d b(t)e 



(58) 



x( 



d^* b d$ d , dt> d d$* b ia 2 ($ d $ b *) 



dRy dR B dR~, dR B 2 dR B dR. 



) 



7 Ollp Ott-y OUQ 



gUJXj 



1 f N 

-N / dnJ2(R0-r k p) 

J k=2 



dU(r k ,-R) 

dRy 



x{Y,cdb{t)e- u ^dn+Y, c Uty tudh ^m- 

d d 

Applying the average procedure i[')l6l) to Eq.(|5"8|). the ob- 
served macroscopic momentum flux (I38[) becomes 



IW (r, t)= J dR/(r - R) £ P* 6 ^ A / dn (59) 



dR^'dR^ + dR^dR^ ~ 2 dR p dR a 



)} 



b ^ fe=2 



dU(r k ,K) 



ORa 



B. Shear viscosity 

Next, we notice that no momentum flux exists in an 
equilibrium state: 

n$(r,i) = Y, P * b f dR f( r H-)n*|(R,t) - 0, (56) 

*6 7 



x{]T c d b(t)e- it ^$ d $ fc * + 2 c db (t)e u ^^* d }. 

d d 

Eq. ([5^|) is first order in H\(t). The order v x0 /l y contri- 
bution comes from Hd- To recover usual viscosity (the 
proportional constant between momentum flux and ve- 
locity gradient v x q/I v ), in Cdb(t), we only need the first 
term of Eq.([i5"|). 
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VI. TIME SCALES: JUSTIFICATION OF THE 
PRESENT SCHEME 

A. Justification of the method of variation of 
constants 

Applying the method of variation of constants (MVC) 
to a macroscopic system worries many scientists 0, [24], 
l57l - [6l"l | . The kinetic approach insists that the long time 
behavior of a large system should be described by a time 
coarse-grained master equation[l, [U [59l - l6l"l | . We will 
show that with some caution, the MVC is applicable to 
calculate shear viscosity. 



1. Wave function description 

We first estimate the time period during which the 
system can be described by a wave function. The thermal 
reservoir plays the role of damping for S. Eq. ([23|) should 
be revised as 



{H + H f + H C + H d )* a - iPM a , (60) 



where P is the net power exerted on S. Usually the 
last term is not written out, but is implicitly taken into 
account S receives the power generated by the exter- 
nal field and velocity gradient, in the same time S delivers 
extra energy to B to maintain the equilibrium with B. If 
at the initial moment t = — oo, the system is exactly in a 
pure state \& Q (— oo) = $ a , *& a (t) is uniquely determined 
by Eq.® and is finite for any t > -oo. Then Eq. (j36)) 
and the corresponding macroscopic equation obtained by 
the averaging procedure (|5l6p works for any t > — oo. 

However, we only know P* a , an incomplete specifica- 
tion of the initial state of <S[3!|, so that a wave function 
description for S cannot last very long. Let us use 



P*°|l 



dt' / dr$:[- 



E 

b&a) 



P q 



dt't 



it LJ a b 



to define a coherence time t c ° h for state a . The expres- 
sion inside the modulus sign on the LHS is the remaining 
probability amplitude of state $ a at time t if \l/ a (t = 
— oo) = the expression inside the modulus sign on 
the RHS is the probability amplitude of state <J> a at time 
t if ty b (t = -oo) = <P b . When t > C^ x = max{t c a oh , Va}, 
it is not possible to describe S with a wave function. The 
energy exchange between S and B makes the system a 
member of the canonical ensemble at temperature T, and 
the system distributes itself into various stationary states 
{$ a } according to Eq.J2]). The incomplete specification 
of initial state (the indeterminacy of initial state) makes 



the system dephase in a time period i^ax- The phase ran- 
domization continues ceaselessly, so that one no longer 
needs the solut ion ^ a (t) of Eq.([23l) for t > t™£ x . 

Expressions (|45|46|) indicate that Dirac's perturbation 
theory can be used only for a time interval much shorter 
than 



(62) 



We can see that t p is the same order as i^ax- I n other 
words, the method of variation of constants fails long 
before the system loses its wave function description. The 
multi-scale method and Fano's ansatzQ can be used to 
obtain ^ a (t) for a longer period: we will discuss them 
elsewhere. As long as the system can be described by a 
wave function (t < f™^), we can apply the method of 
variation of constants to calculate the time evolution of 
wave function. 



Error arising from the indeterminacy of the interaction 
time 



We show that the indeterminacy of the interaction time 
will not lead to any serious error in the results obtained 
by the MVC. 

From 



dt 



i(E b -E a )t/h 



dT^H^a, 



(63) 



we see that the characteristic interaction time is (a;f, a + 
cj) _1 , where oj is the frequency of external conditions. 
For some interaction time t, 



h(uj ba + LU) 



(64) 



Because any external disturbance (including velocity 
(61) gradient) is only exerted on a few degrees of freedom, 
the external distance can quickly relax to other degrees 
of freedom [EH . There is an indeterminacy At in the inter- 
action time, which is the order of microscopic relaxation 
time t m ~ h/U. The error caused by this indeterminacy 



dg, 



dt 



-At 



-ie 



i(E b -E a )t/ 



U 



The relative error is 



« 1. 



(65) 



(66) 



In addition, the phase randomization induced by the en- 
ergy exchange between bath and system leads to a de- 
structive interference in any observable which is bilinear 
in wave function ty a (t), i.e., the error will not accumulate 
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with time. Therefore the indeterminacy in the interac- 
tion time does not cause serious error in any observable. 

This conclusion will not change for two degenerate lev- 
els Eb = E a in a static held ui — and for two levels that 
are in resonance with external condition Uba + uj = 0. 
In both situations, the characteristic interaction time is 
h[J dT$k if i<I> a ] _1 , the relative error of transition ampli- 
tude is If' 1 J dT®lH x & a « 1. For details, see Ref.[5|. 

B. Self-consistency 

We show that the MVC is compatible with spatial 
coarse-graining ([5]). The effect of a monochromatic ex- 
ternal field is simple [13]. An electromagnetic field or 
a gravitational field, affects every molecule at the same 
time. The whole system reaches steady state after a mi- 
croscopic relaxation time h/U. The cutoff wave vector ko 
is determined by the precision of macroscopic measure- 
ment, is not related to any time scale. But for a velocity 
gradient, the length scale of spatial coarse-graining re- 
lates to an intrinsic time scale. 



1. Monochromatic velocity gradient and external field 

If a velocity gradient with time dependence cosut is 
applied to S by a plate or a rod, after a time period sev- 
eral times of the characteristic time t r f = max{(u;ba + 
h[J dT^lHi^a}" 1 }, the molecules within a dis- 
tance t r fh/(md) to the plate or rod reach steady state. 

9ba (*) ano - momentum flux II 7 ^(r, t) contain terms with 
cos ujt factor and terms with smut factor. The coefficient 
of coswt gives us the real part of the viscosity Rer?(w), 
the coefficient of sin cot gives us the imaginary part of the 
viscosity Im?7(u;). For the aim of computing viscosity, we 
choose cutoff wave vector fco through fc^ 1 = t r th/(md). 

For the steady state under a monochromatic velocity 
gradient, the macroscopic response of S at time t corre- 
sponds to the microscopic response at the same time t. 
The averaging procedure (|5I6[) reflects this fact. 

2. Disturbance with arbitrary time dependence 

For a time-dependent mechanical perturbation, the 
macroscopic response at time t depends on all the 
past history (— oo,t) of the external conditions. There- 
fore the macroscopic response is a convolution of all 
monochromatic transport coefficients and the external 
disturbance @. For a velocity gradient with arbitrary 
time-dependence, the macroscopic response is more com- 
plicated. 

For an internal disturbance with arbitrary time depen- 
dence, we make a temporal Fourier resolution: 

/oo 
dujH lu e iut . (67) 
-oo 



Because H\{t) in Eq. (|23|) is complex, we will not require 
H\{ui) = [Hx{—ui)}*. The first order transition amplitude 
is a sum of every monochromatic transition amplitude: 

/OO • ft /> 

du\-~ / dt'e^+^t' / drSlHM 
-OO >/ — OO J 

(68) 

According to Eqs. (|58l38p . the macroscopic momentum 
flux is a frequency integral over all monochromatic mo- 
mentum flux. One must notice that this result originated 
from first order perturbation theory. Unlike a classical 
oscillator, in Eq.(|23|) the external conditions are multi- 
plied by *„(*)• Although Eq.([23J is linear in * Q , the re- 
sponse to Hi(t) is not the sum of the responses to {H\ u } 
in general. 

VII. NON-EQUILIBRIUM DENSITY MATRIX 

A. Equivalence between two definitions of density 
matrix 

The density matrix of the system is often defined as[39| 
1 M 

Pnm = J^J2 a h^ m '^ a (i) fa'*)' ( 69 ) 

where M. is the number of systems in an ensemble, 
d(£)(ra, t) is the probability amplitude that the sys- 
tem is in state <!>„ at time t: 

*(t) = ^o(i,t)*,. (70) 
i 

The density matrix p can also be defined by the initial 
value problem Q: 

ih^ = [H + H 1 (t),p\, p(t=-co) = p e «. (71) 

Here p eq = e~^ H ° /Z is the equilibrium density operator, 
and Z —Tre~^ H ° . By means of 

th^^ = a p {t)E p + Y j a l {t){ P \H 1 {t)\l) ) (72) 

it is easy to check that if H\ it) is hermitian, p is the same 
as p in Eq.flBSJ). 

Integrating by parts, Hd becomes 

H i = -\P*' E ^/dR/tr^-R) 

N 

x ^{Jfo - R) Pjl + Pjl 8(r j - R)}, (73) 

where pj 7 = —ihd/drjy (7 = x,y,z) is the momen- 
tum operator for the j th particle. We see that Hi(t) = 
Hf(t) + H c {t) + Hd is hermitian. Thus we can use either 
Eg. (15^)) or Eg . (fTTj) to define density matrix. The for- 
mer is more convenient for us since we already obtained 
probability amplitudes in Sec IV Al 
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B. Equivalence between procedure (|5I6[) and 
average over density matrix 



The perturbation expansion (|48|) for a time-dependent 
wave function ^ a (t) leads to a similar expansion for the 
density matrix 



Pc b (t)=P^{t)+P^{t) + 



(74) 



where the superscripts indicate the order in H\(t). Be- 
cause (i) the system is enclosed by a bath with tempera- 
ture T; and (ii) a state with initial value $ a is described 
by Ea.(|l8"|). The zero order density matrix isJH, HU 

Pcb(t) ~ P^ b &bc- It is easy to check that p^{t) = 0, 
and the first order correction for off-diagonal element is 



p$(t) = p* b c$(t)+p*'4 1 Ht). 



(75) 



In the basis {$ a }, the matrix element of the momentum 
flux operator (|4"0)) is 



[IW(r)] bc = ($ b |n Q/3 (r)|$ c ) = —N I dn (76) 



d$ c d$* b d<f>* b <9$ c 1 9 2 ($ c $^ 



dr a drp dr a drp 2 dr@dr 



1 f N 
--N J d-rijjfo-r^) 



According to standard statistical mechanics 0, [U, [39| . 
the observed macroscopic momentum flux H a p (r, t) is 

n^(r,i)= /'dR/(r-R)^[n^(R)] bc p«(t). (77) 



Substitute Eq.(|76jl into Ea.([77|). we obtain Eg. ([59]) which 
was obtained by the microscopic response method. 



C. Constancy of the entropy S s of system S 

The entropy of system S is defined by@, HI 

S* 5 ^) = -Trplnp = -^(nlplnpln), (78) 

n 

where (ri • • • rjv|n) = <& n is the eigenfunction of Hq. The 
time dependence of entropy is included in the density 
operator p(t). The time rate of change of S (t) is written 

as 



dS s (t) ^ dplnp 



dt 



Ot 



(79) 



Let A = 1 — p, then p = 1 — A. It is easy to check that 

_dA r ' 



i Ti- 



ck 



[H +H lt A n ], n = 1,2,3, • 



(80) 



Making use of 

\np = ln(l - A) 
Eq. ([80]) leads to a well-known result [3 



„ A 2 A 3 
~2 3~ 



ot 



(81) 



Then, we combine Liouville equation (jTTj) and Eq. (l8"TT) . 

to obtain 



By means of Eq.®, Eq.® becomes 



(82) 



dS s (t) 
dt 



^{n\[H o + H 1 (t),p\D.p]\n) = 0. (83) 



If a system is in good thermal contact with a reservoir 
such that the heat generated by the external field and 
velocity gradient can be instantaneously transferred to 
the surrounding bath, its temporal evolution can be de- 
scribed Eq.(|2"3]l or Eq.([?l"]>. The entropy S s of system 
S is a constant. The irreversibility (entropy produc- 
tion) of the momentum transport process is reflected in 
Eas. d53l52l . 

Eq. ([83l) does not depend on the initial condition for p. 
The key reasons are that (1) the system is described by 
a 'Hamiltonian'. The density matrix satisfies a memory- 
less Liouville equation ([71)1: (2) Tr AB =TtBA is correct 
for two arbitrary operators A and B. For the same rea- 
sons, if the entropy of a microscopic canonical ensemble 
is defined by Eq.fl7S), one has also0| dS s {t)/dt = 0. 
Removing irrelevant degrees of freedom can avoid this 
difficulty!; 



D. Cumulant expansion: connection with previous 
methods 



Eqs. (171l81l82[) have a common form 
3D 

iti— = [H + Hi(t),D], 
ot 



(84) 



where D = p. lnp and p In p. One can solve Eq. ([54l with 
perturbation method^, [2J] 

D{t) = £>(°) (t) + (t) + (t) + ■ ■ • , (85) 

where the superscript indicates the order in Hi (t) . Before 
we introduce external field and velocity gradient, the sys- 
tem is in equilibrium. The initial value of D(t) is its equi- 
librium value D eq : £>(°)(-oo) = D eq and DW(-oo) = 0, 
q = 1, 2, • • • . Since D eq is function of Ho, the zero order 
solution of Eq.flHU isQ 

= e -i(t-t')Ho/h D (0)^y(t-t')H /H = peq (gg) 
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The first order solution isQ 



£>W(t) = ~ / dt'e- t ^ t ' )Ha / h [H 1 {t'),D ec i]e^ t - t ^ Ha ' h . 



It is easy to check that: 

(n\D {1) {t)\n) = 

for any eigenstate $ n of -ffo- 

For D = lnp, [lnp](£ = — oo) = lnp eg , so that 



(87) 



(89) 



[hxp]^{t) = -pH Q -\nZ. 
The first order solution is 

[l np ](i)( t ) = ^ f* dt'e^ t ^ Ha / h [H 1 (t'),H \e^ t ' t '^/ h . 

(90) 

We recognize that [Hi(t'), Hq]/H is the dissipated power 
caused by the external field and velocity gradient. Ex- 
cepting a local non-uniformity of temperature which does 
not appear in the canonical ensemble, [In p] W (t) is Mori's 

-RQ-To obtain the term linear in velocity gradient v x q/I v , 
we only need [Ha, Ho]/h. By means of Eqs. (171251) . the 
term linear in v x q/I v is 

(b\[lnp]W(t)\a) = -P**P f dt> Kit 



J dRfir^-Rme-^-^/* 



2m4- LVr '" v±J " ;j 9r j7 2m^ L " 3 dr., 



[V^fo-R)]- 



]-V, 



'31 



4m 



E[ V ' 3 ^ ( ? R) ]}e^ t,)Ho/h \a). (91) 



Eq.(|9ip is j3 multiplied by the dissipated energy from 
t = — oo to t. 

Eqs. (|89l91l) forms a cumulant expansion for the density 
matrix p = e^ lnp ^: 

p = Z- 1 exp{-f3H Q + pnp] {1) (t) +•••}. (92) 

The operator in {} on RHS of Eq. ([9l"T) is the diver- 
gence of the stress tensor (momentum flux). If we re- 
place the momentum flux operator with its expectation 
value, and change the source of velocity gradient from 
boundary surfaces to bulk, Eq. (|9"2")l is reduced to the non- 
equilibrium density matrix obtained by many different 



methods[25j, |32j, |48|, |49H . We need not invoke Onsager's 
regression assumption [25| and nonconservative force from 
bathJUIli. 

VIII. SUMMARY 

According to the microscopic response method 0, 
I20I [23j . the observed macroscopic momentum density 
Eq. (|13j) in a many-body system can be obtained from 
the microscopic momentum density in a pure state by 
spatial coarse-graining and averaging over all possible ini- 
tial conditions. If we adopt a no-slip boundary condition 
Eq. (ll5[) . we can view a velocity gradient as a constraint 
on the many-body wave function of the system Eq. (1161) . 
The evolution equation for wave function Eq. (|23[) is then 
derived from the Lagrange multiplier method. The La- 
grange multipliers have been obtained by the method of 
variation of constants, Eqs. (j51l52l53[) . They express the 
interactions on the system exerted by the moving bound- 
ary plate and three static plates. 

From the evolution equation Eq. (1231) , we calculated the 
time rate of change of mass density, momentum density 
and total energy in a pure state: Eqs. (134l36l32[) . The dis- 
sipation caused by a velocity gradient contains [P* a ] 2 , 
in contrast to the dissipation caused by a mechanical 
disturbance [2 3] which contains P* a . This is an obviously 
statistical character, velocity gradient is an internal dis- 
turbance. 

Although velocity gradient is an internal disturbance, 
by means of the conservation law of momentum in a 
pure state, we could still define momentum flux in a 
pure state Eq. (|3"7) . The observed macroscopic momen- 
tum flux Eq. (|38[) is obtained by spatial coarse-graining 
the momentum flux in a pure state and averaging over 
all possible initial conditions. The shear viscosity can 
be read out from Eq. (|59")) . Comparing to the traditional 
theories of viscosity, the present ansatz does not need 
temporal coarse graining. Taking a spatial coarse-grained 
average over the microscopic response automatically con- 
tains temporal coarse-graining and coarse-graining in the 
eigenvalue spectrum of collective variables. We com- 
pared the Hamiltonian derived from Lagrange multi- 
plier method with the phenomenological non-equilibrium 
Hamiltonian. The non-equilibrium density matrix im- 
plied by the present method can be reduced to those 
derived by other methods. 

The present ansatz can be generalized to all internal 
disturbances. We can view concentration gradient, tem- 
perature gradient and velocity gradient as constraints on 
the many-body wave function of system, and have an uni- 
fied theory for diffusion, thermal conductivity and viscos- 
ity. This work is in progress. 
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